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Abstract. We are interested in the life span and the asymptotic behaviour of the solutions to a 
system governing the motion of a pressureless gas, submitted to a strong, inhomogeneous magnetic 
field e~^B(x), of variable amplitude but fixed direction — this is a first step in the direction of the 
study of rotating Euler equations. This leads to the study of a multi-dimensional Burgers type system 
on the velocity field u^, penalized by a rotating term e~^Ue f\B{x). We prove that the unique, smooth 
solution of this Burgers system exists on a uniform time interval [0, T]. We also prove that the phase 
of oscillation of is an order one perturbation of the phase obtained in the case of a pure rotation 
(with no nonlinear transport term), e^^B{x)t. Finally going back to the pressureless gas system, we 
obtain the asymptotics of the density as e goes to zero. 

Resultats asymptotiques pour la magneto-hydrodynamique sans pression 

Resume. On s'interesse au temps d'existence et au comportement asymptotique des solutions d'un 
systeme modelisant le mouvement d'un gaz sans pression, soumis a un fort champ magnetique £~^B{x), 
d'intensite variable mais de direction fixe — cela etant un premier pas dans la comprehension des 
equations d'Euler en rotation rapide. Cela conduit a I'etude d'un systeme de type Burgers multi- 
dimensionnel sur le champ de vitesse u^, penalise par un terme de rotation e~^u^ f\B{x). On demontre 
que la solution unique reguliere de ce systeme de Burgers existe sur un intervalle de temps uni- 
forme [0,T]. On montre aussi que la phase d'oscillation de est une perturbation au premier ordre 
de la phase obtenue dans le cas d'une rotation pure (sans terme de transport non lineaire), e^^B(x)t. 
Enfin en revenant au systeme des gaz sans pression, on obtient le comportement asymptotique de la 
densite quand e tend vers zero. 

1. Introduction 

The aim of this paper is to study the asymptotic behaviour of a fluid submitted to a strong 
external inhomogeneous magnetic field. 

The case when the field is constant has been studied by a number of authors, both for com- 
pressible and incompressible models of fluids (see for instance Q , or [7] for incompressible 
fluids, and 4, or JU] for rarefied plasmas). In that case, one can not only derive the asymp- 
totic average motion (which is given by the weak limit of the velocity field), but one can also 
describe all the oscillations in the system and possibly their coupling : the filtering techniques 
used for that rely on explicit computations in Fourier space. 
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In the case when the magnetic field is inhomogeneous, those methods are not relevant any 
more. Weak compactness and compensated compactness arguments allow nevertheless to 
determine the average motion (see [B] in the case of a rarefied plasma governed by the Vlasov- 
Poisson system, and in the case of a viscous incompressible fluid). In order to describe 
the oscillating component of the motion, one has to understand the interaction between the 
penalization and the nonlinear term of transport : indeed one expects that the flow modifies 
substantially the phase of oscillation (which is of course inhomogeneous). 

We propose here to analyse this interaction for a simplified model of magneto- hydro-dynamics, 
the so-called Euler system of pressureless gas dynamics. 

1.1. A simple model for magneto- hydro-dynamics. We consider the following system 
of partial differential equations : 

dtp + V- {pu) = 0, xeK^,t>0 

(1.1) dtipu) + V- {pu ^ u) = pu A B , xGR^,t>0 

p{t = 0) = po, u{t = 0) = no X G R^, 

where p denotes the density of the fluid, u its mean velocity and B is the external magnetic 
field (V- B = 0). The first equation expresses the local conservation of mass, while the 
second one gives the local conservation of momentum provided that there is no internal force 
(no pressure). This assumption is relevant only in some particular regimes (corresponding 
to sticky particles j2j). Prom a physical point of view, this may seem a strong restriction, 
but it allows to perform a first mathematical study of that type of inhomogeneous singular 
perturbation problem : indeed in this special case a major simplification arises since the 
equation on the mean velocity can be (at least formally) decoupled from the rest of the 
system : 

dtu+ {u -Vju = u A B xGR^,t>0. 

We then obtain a system of Burgers' type, that is a prototype of hyperbolic system. A work 
in progress should extend the present results to more realistic models, in particular to the 3D 
incompressible Euler system. 

In order to further simplify the analysis, we assume that the direction of the field B is constant 
B{x) = -b{xi,X2)e3, (xi,X2)gR^ and 63 = *(0, 0, 1) 

which allows to get rid of the geometry of the field lines (for detailed comments on this 
subject see for instance j^. Remark 1.4). Any solution to the system (jl.l|) has then uniform 
regularity with respect to the variable X3. To isolate the phenomenon of inhomogeneous 
oscillations with instantaneous loss of regularity, we restrict therefore our attention to the 
2D singular perturbation problem, in the plane orthogonal to the magnetic field. We finally 
have : 

dtp + V- {pu) = 0, xGR^f>0 

(1.2) dtipu) + {pu(S)u) = ^pu^, xGR^t>0 

p{t = 0) = Po, u{t = 0) = no, X G R^, 
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where u denotes the vector field with components (u2, —ui), and the intensity b of the 
magnetic field satisfies the fohowing assumptions : 

b e C7~(r2) n vf2'°°(r2), {m) 

inf b{x) = 6„ > 0. {HI) 



A standard fixed point argument then allows to prove the local well-posedness of (|1.2|) . The 
result is the following. 

Theorem 1. Consider a function b satisfying assumptions {H0){H1). Let po be a nonnegative 
function and uq be a vector-field in H^(R?) (s > 2). Then, for all e > 0, there exist S 
]0,+oo] and a unique solution of {pe,Ue) G L'^^{[0,Ts[, H'{R'^)). 

Note that the lifespan of the solution depends on e, and that the lower bound on coming 
from the Duhamel formula goes to zero as e ^ 0. The first difficulty to study the asymptotics 
e — > consists then in understanding why the magnetic penalization does not destabilize the 
system, and in proving that the solution {p^,Ue) exists on a uniform interval of time. 

1.2. Formal analysis. Before stating more precise results on the lifespan of the solutions 
and on the asymptotics e — > 0, we have chosen to give some simple observations about the 
problem to guide intuition. In this first approach we restrict our attention to the analysis of 
the equation governing the velocity. 



The first step of the formal analysis consists in determining the mean behaviour of the velocity 
field, that is the weak limit of u^. We have 

Ms = ^ {dtUe + {Ue ■ V)ne)"^ . 

As b is bounded from below, if we are able to establish convenient a priori bounds on u^, this 
will imply 

in some weak sense. This means that we expect the velocity to oscillate at high frequency (on 
vanishing temporal or spatial scales). 



Another way to get an idea of the asymptotic behaviour of the velocity is to study the simple 
case when b is constant. The group of oscillations generated by the magnetic penalization is 
then homogeneous : 

fbt\ ^ . fbt\ 
K \- j u = u cos I — ) ~ ^ I — ) ' 

which corresponds to the rotation with frequency 27r6/e. As the coefficients are constant, this 
group is not perturbed by the transport. Classical filtering methods (see namely IZIiIH]) can 
then be applied : setting 

= R{ — ] Ue 



leads to 

dtVe + Q { -,Ve 
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where Q ., .) is a quadratic form with bounded coefficients depending on t/e. As there is 
only one oscillation frequency, there is no resonance, which implies that 

Ve Uq 

in some strong sense, provided that convenient a priori bounds on (and consequently on v^) 
hold. This means that we can describe completely the oscillations and get a strong convergence 
result. Of course, we get as a corollary that the lifespan is uniformly bounded from below, 
and we even expect that +00 as e — > 0. 



The case we consider here is much more complicated. The group of oscillations generated by 
the magnetic penalization is again very easy to describe : 

A \ fb{x)t\ ^ . fb{x)t 
K \-,x j u = u cos I / ~ ^ ^'^^ \ 

but it is non homogeneous, which entails 

• a loss of regularity (R u blows up in all Sobolev norms H^CR?) for s > 0); 

• an interaction with the transport operator (with the same definition of = R[—^,x^ Ug as 
previously, we do not expect dtv^ to be bounded in any space of distributions). 

The stake behind this model problem is to understand how to overcome these difficulties. The 
first step is to explain how the phase of oscillations is modified by the fiow : note that even a 
small correction on the phase changes strongly the vector field. Then we have to establish a 
strong convergence result using a new method : classical energy methods fail because of the 
lack of regularity on approximate solutions. Here an appropriate rewriting of the system by 
means of characteristics associated with the flow allows to understand the underlying structure 
and to answer both questions : in particular we will see that the spaces which are well adapted 
for this type of study are constructed on L°°(R?). In the case of incompressible dynamics, the 
analysis will be therefore much more difficult since the transport is replaced by a non-local 
pseudo-differential operator. 

1.3. Main results. As long as the solution {pg,Ue) of system ()1.2|) is smooth, the velocity 
satisfies the equation of Burgers type 

dtUe + {Ue ■ V)Ue + -U^ = 0, X G R^, t > 

(l.oj e 

Ug{t = 0) = X G 

Using refined a priori estimates on this last equation, we can prove that for all e > it admits 
a smooth solution on a uniform time T > 0. We will prove the following result. 

Theorem 2. Consider a function b satisfying assumptions (HO) (HI). Let pQ be a non- 
negative function in P^*~^'°°(R^), and uq be a vector-field in W^'°°{'R'^) (s > 1). Then, 
there exists T* g]0, -l-oo] such that, for all j" < j"* gjjc! all £ < et, there is a unique 
{ps,Ue) G L~([0,T], VF^-1'~(R2) X VF^'°°(R2)) solution of (which is nevertheless not 
uniformly bounded in L°°([0, T], VF^-1'°°(R2) x Ty^'°°(R2)) for s > 0). 

In this framework, it is relevant to consider the asymptotics e ^ on the time interval [0, T]. 
The same type of computations as previously allows to prove that the velocity field behaves 
almost as in the constant case (with slight modifications of the phase of oscillations). 
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Theorem 3. Consider a function b satisfying assumptions (HO) {HI). Let uq he a vector- 
field in W^''^{R'^) (s>l). For all T < T* as in Theorem\^and all e < et, denote by Us the 
solution of C3) in L°°([0,r], VF"'°°(R2)). Then, 

Ue{t, x) — (^q{x) cos 9e{t, x) — Uq {x) sin ^^(t, x)^ 

converges strongly to in L°°([0,T] x R^), where the phase 9^ is defined by the following 
equation 

b(x)t 

(1.4) 0^{t, x) = -y^ - tuo{x) ■ V log b{x) sin 6^(1, x) 

+tu^{x) ■ V log b{x) cos x). 

Rewriting the equation on the density ps with a transport term and a penalization term 
(coming from the divergence of which is of order 1 /e) 

, ^ dtPe + ■ ^Pe + Pe^- Ue = X G R^, f > 

(1.5) „ „ 

Pe(t = 0) = p°, x€R2, 

we can then determine the global asymptotics of the Euler system of pressureless gases H1.2() . 

Theorem 4. Consider a function b satisfying assumptions {H0){H1). Let po be a nonnegative 
function in W'^'°°{B?), and uq be a vector-field in W'''°°(R'^) (s > 1). For all T <T* as in 
Theorem 12 and aU e < et, denote by (p^, Ue) the solution of in L~([0, T], W'-^^°^{Y{?) x 
Ty^'~(R2)). Then, 

Pe{t, x) — poix) ^1 + tuQ ■ V log(x) cos 6e{t, x) — tuQ ■ V log b{x) sin 6eit, x)^ 

converges strongly to in L°°([0,T] xR^), where the phase 9^ is defined as previously by i|I.4|) . 



Let us comment a little on the proof of those theorems, and give the structure of the paper. 

It is quite clear that energy methods will not enable us to have a good control on the asymp- 
totics of (pejUs), since as soon as we want a control on derivatives of unbounded terms 
will appear. So the most appropriate way to study System (|1.2|) is to rewrite it using the 
charcateristics of the flow and to study those characteristics precisely. 

Section 2 is therefore devoted to rewriting System ()1.2|) in characteristic form, and in the 
derivation of a few a priori estimates. 

In order to establish the existence of a solution (/>£,-Ue) to system p.2() on a uniform time 
interval [0, T], it is enough to see that the solution is well-defined (and smooth) as long as the 
flow generates a diffeomorphism ^^(t, .) 

^^(t,x) = Us{t,Xs{t,x)) 

hence to prove that the characteristics cannot cross before time T. The precise estimates 
on DXe leading to Theorem 2 are performed in Section 3, they use in a crucial way some 
results of non-stationary phase type. 

The asymptotic behaviour of Us{t, X£{t, .)) and p^{t, X^{t, .)) is then simply obtained from 
the explicit approximation of the characteristics X^, using Taylor expansions for the various 
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fields. In order to establish the convergence results stated in Theorems |31 and 01 the main 
difficulty is therefore to get a precise description of the inverse characteristics X~^(t, .), which 
is done in Section 4. 



2. Appropriate formulation of the system 

As pointed out in the introduction, energy estimates do not seem to be the right angle of 
attack for our problem. We shall therefore in this short section present a new formulation of 
System (|1.2|) . by means of characteristics (Paragraph In that way some a priori estimates 
can be deduced immediately (see Paragraph 12. 2|) . 



To simplify notation, from now on we shall drop the index e in and simply write u (and 
similarly for any other e-dependent function). 



2.1. Trajectories associated with the flow. Let us write System (|1.2|) in the following 
form: 

dX 

Ik 



u{t,X), X\t^o = x 



(2.1) ^(p{t,X))+pV-u{t,X)=0, p\t=o = Po 

at 

±(uit,X)) + ^u^{t,X)=0, n|i=o = no. 
at e 

As seen in Theorem^ there is a solution to System (|1.2() for a time depending on e, and as 
long as the trajectories do not intersect we can write in particular 



(2.2) u{t, X{t, x)) = uo{x) cos _ u^(x) sin 



where we have defined the functions 

^{t,x)=[ P{s,x)ds, P{t,x)=h{X{t,x)) 
Jo 

(these functions are well defined as long as the characteristics do not cross each other). 

If u is smooth enough, then p is uniquely defined by the transport equation it satisfies. So 
from now on we can concentrate on u (and X). As one of the aims of this article is to prove 
Theorem 1^ (which will be achieved in the next section), we shall from now on call the 
largest time before which no characterestic intersect; one of our goals is to prove that is 
uniformly bounded from below as e goes to zero. 

In the next paragraph we are going to derive from (|2.1|) and (|2.2jl some easy a priori estimates 
for times < t < T^, which will help us prove Theorem [21 in the following Section (31 and 
Theorems O and ^ in Section 



2.2. A priori estimates. Formula (|2.2|) immediately enables us to deduce the following a 
priori estimate: 

(2-3) lkllL°°([o,T'=[xR2) < 2||mo||l°° , 



which imphes that 
(2.4) 
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dX 



dt 



< 2||no IIloo . 

L°°([0,T=[xR2) 

In particular X — x remains bounded in space for all times < t < T"^, and we have 
(2.5) Vt G [0,r^[, ||X(t,-) -x|lioo(R2) < 2t||no||ioo(R2). 

Since (3{t,x) = b{X{t,x)), we have 

\\dtP\\L°°iio,T'^[xn?) < I|V6||loo(r2) 

(2.6) 



dX 



dt 



L°°([0,T=[xR2) 

^ 2||V6||2,oo(p_2)|[no||L°°, 
as well as 

(2.7) VtG[0,T^[, VxGR^ 6„ < /5(t, x) < ||&||ioo(R2) 

with 6_ defined in (HI). 

Now we are going to look for an approximation of X: integrating formula (|2.2() in time yields 



(2i 



X{t,x) = X + uo{x) J cos ^ ^ — J ds — uj^(x) J sin 



(s,a;) 



recalling that (pit^x) = / 6(X(s,2;)) ds. The following section will be devoted to a precise 

^0 

study of the trajectories X, which will enable us to infer Theorem |2l 



3. Study of the trajectories 

Formulation 1)2.11) of System (|1.3|1 shows that the study of the Euler system of pressureless 
gases with magnetic penalization comes down to a precise analysis of the characteristics, and 
in particular of their invertibility. 

In this section we will establish that the trajectories defined by 1)2. 8() are invertible on a time 
interval [0,r=[ with 

lim T" = T* > 0, 

e^O 

where T* depends on the magnetic field b and on the initial velocity field uq. This result is 
based on an asymptotic expansion of the Jacobian 



def 



which implies that 



J{t,x) = \det{DX{t,x)) I, 



yte[0,T*[, liminf J(t,x) > 0. 

£^0 



The asymptotic expansions of X and DX (Paragraphs 3.2 and 3.4) are obtained using some 
results of non-stationary phase type and the L°°-bounds established in Paragraphs 3.1 and 3.3. 
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3.1. Bounds on X(t,-). The first step of the analysis consists in showing that for any 
point X G R^, the characteristic stemming from x stays in a ball of size 0{e) around x. 
This shows that the rotation has a drastic influence over the transport by u. 

We have the following proposition. 

Proposition 1. Let x G be given, and let X{-,x) he the trajectory starting from x at 
time 0, defined by i|2.8l) . As long as it is defined, it satisfies 



Mt < min(r, 



\X(t,x) -x\< 4— ||uo||l°o(1 + T 



Proof of Proposition ^ The proof is an immediate application of the non-stationary 
phase theorem. As we will be using such arguments many times in the following, let us state 
and prove the following lemma, which will be invoked systematically in the next sections. 

Lemma 1. Let T he a given real numher, possihly depending on e. Let F he a function 
uniformly bounded in VF^'°°([0, T], L°°(R^)), and let j3 he a positive function, also uniformly 
hounded in W^'°°{[0,T], L°°Ir'^)), and hounded hy helow hy 6_. Then for all t € [0,r] and 
all X £ R^, the following bounds hold: 



and 



F(s, x) cos 



Pis',x) 



F{s, x) sin 



(3{s',x) 



ds I ds 



< £ 



ds I ds 



lint, Oil 



L°°(R2) 



+ t 



(3{s, 



/ ||F(t,-)||^^(R2) + ||F(0,.)||^^(R2) ^ ^ 



d. 



Pis, 



L°°([0,t]xR2)^ 



L°°([0,t]xR2)^ 



Proof of Lemma ^ The proof is a simple application of the nonstationary phase theorem: 
an integration by parts leads to 

t y fJ{hX) \Jo 



F(s, x) cos 







and similarly 

F(s, x) sin 



I3{s',x) 







Pis',x) 



(3{t,x 
Jo 



ds 



e 



ds' ] ds = e 



F{t,x) 
e—, cos 



F{^,x) 

/3(0,x) ^ I3{t,x) 
, F{s,x) 



I3{s',x) 



ds' I ds, 



P{s,x) 
e 

e 



ds 



ds' I ds. 



The result follows immediately. 
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Now let us go back to the proof of Proposition E Recalling formula (|2.8j) , we simply apply 
Lemma^to the case F(t,x) = uo{x) to get 

dsl3{s, ■ 



\X{t, x)-x\< 4£^ML^ + 2st 



Uq- 



/32(s, 



L°°([0,t]xR2) 



Estimates (|'2.6() and 1)2 .Tf) immediately yield Proposition ^ 



3.2. Asymptotics of X{t, •). The same type of computations based on the non-stationary 
phase theorem allows actually to obtain an explicit approximation of the characteristic X at 
any order with respect to e (in fact we will stop at order 2 but the argument can be pushed 
as far as wanted if necessary). 

Lemma 2. For any point x £ and any time t < min(r, T^), the following approximation 
of the trajectories defined in i|2.8|) holds: 



X{t, x) — X — e 



uojx ) 
h{x) 



■ sm 



^)+e^ (l-cos(^) 



+ £tv{x) 



< Cte' 



where the drift velocity is given hy 
and Ct denotes a constant depending only on T, uq and b. 



(n^-V6)uo(x)-(uo-V6)4(x) , 



Proof of Lemma [2J Let us write the following expression for X{t,x), obtained from 1)2. 8|) : 

we have 

X{t,x) = x + R%t,x), 

with 



(3.1) 



W{t,x) 



def 



uq{x) I cos 



(s,x) 



ds — 1*0" (2;) 



sm 



{s,x] 



ds 



Rl{t,x) + Rl{t,x). 



We shall only compute the approximation for Rf{t,x), and we leave R2{t,x) to the reader. 
By an integration by parts we have 



(3.2) 



Rl{t,x) = e^sin^^^'^''^ 



£Uo{x) / a 



1 



l3{s,x] 



I3{t,x) \ e 

The first term is easy to approximate : we have, due to Proposition^ 
1 1 



sm 



(s, x) 



ds. 



(3.3) 



Pit,x) b{x) 



< \\Vb\\Loo\X{t,x)-x\ 

4||V6||L°°||no||L°° ||V&||L°°lko||L°° x 
- bl ^' 



for all t < min(T, T*") and all x £ R^. So (3{t,x) can be replaced by b{x) in the first term 
of Rf in ()3.2|) . up to a remainder eTZ^ with ||'7^^||loo([o,t]xR2) ^ Ct- 

Now we need to approximate the second term. Using the fact that 

dsp{s,x) = {u-Vb){s,X{s,x)) 
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with 



ti(s, X{s, x)) = Uo{x) COS 



Uq (x) sin 



(s,x) 



we can therefore write 



5, 



1 



P{s,x) 



sm 



Note that similar computations lead to the following formula, which is useful to estimate iZf- 



\7b(X(s,x)) (2(i)(s,x) 



ds 



/3(s,x) 



cos 



(s,a;) 



ds 



* ^ V6(X(g,x)) . (24>{s,x) 
Uq [x) ■ sm ' 



(3.5) 



Jo 



Vb{X{s,x)) 
262(X(s,x)) 



1 + cos 



262(X(s,x)) 
2(/)(s, X 



ds 



ds. 



Both formulas H3.4() and H3.5() show that new harmonics have been created by the coupling in 
the equation. 

Let us go back to the estimate of the right-hand side in p. 4(1 . To estimate the oscillating 
terms, we use Lemma Q with 

^/ N /X Vb(Xis,x)) , , , , , V6(X(s,x)) 



262(X(s,x)) 



262(X(s,x)) 



We get 



/ Fi(s,x)sin 
Jo 



2(t){s,x) 



ds 



< 2e||uo||L= 



IV6I 



hi 

+et\\uo\\L°^ 



5, 



Vb{X{s,-)) 



b^Xis,.)) 



L°°([0,t]xR2) 



and similarly 

r2(S, X) cos 







e 



ds 



< e\\uo\\L°° 



bl 

-|-et||tio||L°o 



a, 



Vb{X{s,-)) 



L°°([0,t]xR2) 



By and (|^ we have 
Vb{X{s,-)) 



5, 



b'iX{s,-)) 



L°°([0,t]xR2) 



<2||D26|U.™^+4||V6||i.^ 



Plugging that estimate along with (|3.3j) into the definition of Rf in (|3.2|) . we get finally 



i?f (t, x) = 6 ^°/^'^ sin 
6(x) 



f3{s,x) 



ds +e^n''{t,x) - euo{x) 



* u^ix) ■ Vb{X{s,x)) 
2b^{X{s,x)) 



ds. 
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^b{X{s,x)) , Vb{x) . , ^ 

JNow we can approximate ,„ , „ . by , . , up to a remainder e/Z . bo we nave 

o^(A(s,xj) b[x) 

, , /■* Uff(x) • V6(X(s,x)) , u):(x)-\Ib{x) , , 

The estimate of B^{t^x^ is similar and left to the reader. This ends the proof of Lemma|2j 

3.3. A priori estimates on DX{t, •). A necessary and sufficient condition for X(t, •) to be 
invertible is that 

J{t,x) = \det{DX{t,x))\ 
does not cancel. In order to obtain a lower bound on the time (before which the char- 
acteristics do no cross each other), we therefore need to study the behaviour of the deriva- 
tives DX{t, •). First of all we derive a uniform L°°-bound which will allow to neglect some 
terms in the asymptotic expansion. 

Lemma 3. Let x G be given, and let X{-, x) be the trajectory starting from x at time 0, 
defined by As long as it is defined, it satisfies 

yt < min(r,T^), Vx G R^, \\DX{t,x)\\ < Ct, 

where Ct denotes a constant depending only on b, uq and T. 

Proof of Lemma |3J Differentiating 1)2. 8() . leads to 

DX{t, x)-Id = Duo{x) cos (^^iEl^^ ds - Du^ix) sin (^^^El^^ ds 

— uq[x) / L>(p(s, X) sm \ 1 ds Uq [x) / L>(p[s,x) cosy 1 as, 



with 



D(l){s,x) =^ / DX{t,x) ■ Vb{X{T,x))dT. 
Jo 



Applying the Fubini theorem to both last terms, we can set this identity in a suitable form 
to get a Gronwall estimate 

DX{t, x)-Id= Duo{x) cos (^^i^^ ds - Du^ix) sin (^^i^l^^ ds 

(3.6) -^uo{x) J\dX{t, x) ■ V)6(A(t, x)) J' sin J^^i^^ ds dr 

-^u^{x) J\dX{t, x) ■ V)b{X{T, x)) J' cos (^^^El^^ ds dr. 

From formula (|2.8|) we deduce that 

uo{x) [\m(^^^) ds + u^{x) j\os(^^^^^ ds = {X{t,x)-X{T,x))^. 



^ e 

Plugging this identity back into 1)3. 6|) leads to 



DX{t,x)-Id= Duoix) [ cosf^^) ds-D4{x) [ sin f 
(3.7) 1 /■* Jo \ ^ J 

— / {X{t,x) - X{t,x))^ (g) iDX{T,x) ■V)b{X{T,x)) dr 
^ Jo 



ds 
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As in the proof of Proposition^ Lemmanyields the following estimate : for all t < min(T, T^), 



(3.8) 
with 



{s,x) 



ds 



< CteW^UqWlcx, , 



o_ V b- 



From (|3.7j> we then deduce an inequality of Gronwall type 
(3.9) 

\\DXit,-)\\Loo <l + CTe\\Vuo\\L-o + [ \\DX{T,-)\\Loo\\Vb\ 

Jo 

By Proposition ^ 

1 



L°° 



{X{T,x)-X{t,x)) 



dT 



yt < min(r,T''), 



-{X{t,-)-X{T,-)) 



< 2Ct\\uo\\l'x', 



loo 



hence 



\\DX{t, OIIl- < (1 + CTe||Vno||Loo) exp (2C7t|| V6||loo ||Mo||L-i) , 
which is the expected estimate, proving Lemma |21 



3.4. Asymptotics of DX. In view of the results established in Lemma[21 we expect actually 
the derivatives diX{t, x) to behave asymptotically as 

A(t, x) + fi{t, x) cos [ ^^-^ — ^ ] + i/{t, x) sin ' '^^ ' ^ 



where A, /j, and u denote some functions which do not depend on e. Such an asymptotics can 
be justified using the same techniques as in the previous paragraph: let us prove the following 
lemma. 

Lemma 4. Let x € be given, and let X{-,x) be the trajectory starting from x at time 0, 
defined by i)2.8|) . Tfien, for all t < mm{T^ ,T) and for all x G R^, 



DX{t,x) - Id -tuo®V logb cos ( ^i^li^ j + tug (g) V log 6 sin ( ^^'^'^^ 



<CTe, 



where Ct denotes a constant depending only on b, uq and T. 

Proof of Lemma 131 Denote by g the function defined on [0,T] x R^ by 

g(t,x) DX{t,x) - Id - tuo ®V \ogbcos ^^i^^ (g, viogfesin (^^^^ 

In view of (|3.9() . we expect g to satisfy a Gronwall inequality of the following type 

1 



e 



-(X(t,.)-X(r,.)) 



\\Vb\\L^dT + CTe, 



(3.10) \\g{t,x)\\ < / \\g{T,x) 

Jo 

for all t < min(T^,T), where Ct denotes a constant depending only on T, uq and b. 
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Let us postpone the proof of this inequahty for a while, and show how it enables us to infer 
Lemma 13 It is easy to see that 

g{0,x)=0. 

Applying the Gronwall lemma and using Proposition ^ as in 1)3. 9p leads to 

V^<min(^^^), VxGR^ \\g{t,x)\\ <CTe. 



Now let us go back to the proof of (|3.9|) . We first compute 
def 



(3.11) A{t,x) 
By Lemma 121 



-{x{t,x)-x{T,x)y 



{X{t, x) - X{t, x)) O {g{T, x) ■ Vb{X{T, x))) dr. 



— ^(xj sm 

V V e . 

--(x)fcos^ 



sm 



cos 




(n^ • yh)u^{x) + {uq ■ Vb)uoix) I + eTZ^t, r, x), 



t - T 
~W{x) 

where TZ'^ is uniformly bounded in L°°([0,T]^ x R^). Plugging this formula back into the 
integral (|3.11|) leads to 



(3.12) 
with 

and 



-{X{t, x) - X{t, x)) ) (g) {g{T, x) ■ Vb{X{T, x))) dr = Ai{t, x) - A2(t, x), 





def 



1 



Ai{t,x) = J \^-{X{t,x) - X{t,x))] ®{DX{t,x) ■Vb){X{T,x))dT 



A2{t,X 



def 



-p-(xi I sm 
6 \ \ e 



sm 



--(x)(cos(^ 



cos 



^-j-(t - r) [{u^ ■ V6)n^(x) + (uq • V6)no(x)) + e7^"(^, r, x) 



262 (x 

V6(X(r, x)) + (no(x) • V6(X(r, x)))V log 6(x)t cos 



r, X 



(uo (x) • V6(X(r, x)))V log 6(x)t sin 



r,x 



dr. 



From (|3.7() and (|3.8() we deduce that 

^i(t,x) = Id - DX{t,x) 

up to terms of order e. 

In order to estimate the second term A2{t, x), we use again a non-stationnary phase theorem. 
Since the trajectories lie in balls of size e, 

|V6(X(t,x)) - V6(x)| < CrWD^bllL^e, 
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for all X G and all r < t < min(T'',T). Then, as dgfi is uniformly bounded according 
to (|2.6j) . Lemma n shows that 



A2{t,X 

ft /„,± 

+ 

ft 

+ 



I'M 
IM 



t-^(x) sm 



(j){t,x) 



uo ( (t){t,x] 

t—(x)cos 

b \ £ 



Vb{x) 



(x) sin 

, yx) cos 
V " 

Using the identities 



£ 
£ 



{tu^{x) ■ V6(x))V log 6(a;) sin 
(rno(x) • V6(x))V log b{x) cos 



\T,X) 
£ 

(t,x) 



dr 

dT + £Vf{t,x) 



cos^ = 1(1 + cos(2(/>)), sin^ = 1(1 - cos(2(/))), 



we then obtain that 

A2{t,x) = 

up to terms of order e. 
Then (|H.12j) can be rewritten 
1 



( (l){t,x) 

t^(x) sm 

6 V e 







Vb{x) 



-(X(t, x) - X{t, x)) {g{T, x) ■ Vb{X{T, x))) dr = -g{t, x) + £71" {t, x), 



which implies immediately (|3.1U|) and yields Lemma |3 as explained above. 

3.5. Existence on a uniform time intervaL As an immediate corollary of Lemma ^ we 
obtain that X{t, •) is a diffeomophism of on a uniform time interval. Indeed, DX is 
invertible as long as 

\\DX -Id\\Lo- < 1. 

Corollary 1. Consider a function b satisfying assumptions {H0){H1). Let {po,uq) be re- 
spectively a nonnegative function of VF*~^'°°(R^) and a vector-field of Ty*'°°(R^) (s > 1). 
Then, for allT < WuqW^L \ \^b\\iL, there exists > such that System admits a unique 
solution {pe,Ue) G L°°([0,T], iy*-i'°°(R2) x VF^'°°(R2)) for aU e < £t. 

Proof of Corollary ^ By Lemma [IJ the trajectories defined by 1)2. 8|) are continuously 
differentiable and satisfy for all t < T and all x € R^, 



DX{t, x) — Id — tuQ ® V log 6 cos 



(hit, x)\ I „ , , . / (h{t, x) 
' ' ^ + tu^ (g) Vlogfesm ' ' ^ 



< Ct£. 



def 



This implies in particular the following estimate on the Jacobian J{t,x) = \ det{DX{t,x))\ : 



J{t, x) — 1 — tuQ ■ V log b cos 



(l){t,x) 



tuQ ■ V log b sin 



< Ct£. 



Then for T < \\uo\\i^L\\Vb\\^L, there exists £t such that 

ye<£T, yt£[0,T], sup |J(t,x) - 1| < 1, 

which means that X is a C^-diffeomorphism of R^ 
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Moreover, from formula (|2.8)) we can deduce by induction that X{t, .) (and consequently its 
inverse X~^{t,-)) is smooth, its regularity being the same as the regularity of the initial 
velocity field uq. Then the vector field u given by 



belongs to L°°([0, T], W^'°°{'R?)) and it is easy to check that it satisfies System in strong 



The density p is then obtained as the strong solution of the linear transport equation 

dtp + u ■ Vp + pV ■ u = 

whose coefficients belong to L°°([0, T], iy'*-i'°°(R2)), with initial data in Ty^-i'°^(R2). It 
therefore stays in L°°([0, T], H^''-1'°°(R2)). We emphasize once again that no uniform bound 
on {p,u) is available in L°^{[0,T],W'^'^(R'^) x VF"'°^(R2)). 

Theorem 121 is proved. 

Remark 1. The supremuni of the life span of the solutions corresponds to a crossing phenom- 
enon, to be compared with the caustic in geometrical optics. Beyond this time, the differential 
system 



with initial data (x, uo(3^))xgR2 ^^^^^ admits a unique smooth solution, but the application 
(X{t,x),(^(t,x)) X(t,x) is no longer injective, it cannot be lifted. The hyperbolic sys- 
tem no longer has a solution. 



Let T < T* = ll'UollZ-llV^llLi' be fixed. Then, for any e < £t as in Corollary ^ the 
solution {p,u) of System (dHJ with initial data ipo,uo) G W^-^'°°{R'^) x VF'''°°(R2) belongs 
to L°°([0,r],VF*-i'°°(R2) X M/"'°^(R2)). Then it makes sense to study their asymptotic 
behaviour as e — > 0, and the aim of this section is to prove Theorems 01 and 01 

Paragraph 14.11 is devoted to the asymptotics of u{t, X{t,x)) ad p{t, X{t, x)). The last para- 
graph consists in inverting the characteristics in order to infer Theorems and 01 

4.1. Asymptotics of u{t, X) and p(t, X). From the characteristic formulation of System (|1.2|1 
and the asymptotic expansion of X{t, •) we immediately deduce the asymptotic behaviour of 
u{t,X{t, •)) and p{t,X{t, ■)). 

Proposition 2. Consider a function b satisfying assumptions {H0){H1). Let uq be a vector- 
field in iy'°°(R^) (s > 1). For all T < T* and e < st as in Tiieorem|2 denote by u the 
solution of 113) in L°°([0, T], T^^'°°(R2)). Then 




sense. 




4. Study of the asymptotics of and p^ 



u{t,X{t,x)) - [Uo{x)cos{cl)e{t,x)) - Uq (x) sin((/)e(f, x)) 
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>2 



converges strongly to in L°^{[0,T] x R ), at speed 0(e), where the phase (ps is deEned by 
(4.1) 



Mt, x) = ^ - t (u^{x) ■ V) log b{x) 



£ 

Proof of Proposition [2J Let us first recall that 



u{t,X{t,x)) = Uq{x) I COS 

where the phase (p is given by 







e 



ds — Uq{x) / sin | - 



{s,x) 



ds, 



(/)(t,x) = I b{X{s,x))ds. 
Jo 



Then in order to establish Proposition [21 we have to approximate the phase. By Lemma |21 



b{X{t,x)) =b{x) + 



£Uo{x) . ( 4)(t, x) 



■ sm 



e 

eu^{x) 



Vb{x) 



b{x) 



1 — cos 



(t){t,x) 



Vb{x)+e^TV{t,x 



noticing that v ■ Vb = 0. It follows that 
u{t, X(t, x)) = uo{x) cos 



b{x)t ^ uo{x) ■ Vb{x) 



b{x) 

u^{x)-Vb{x) /■* / 



b{x) 



1^ sin ds 

[s,x) 



-Uq{x) sin 



b{x)t uq{x) ■ Vb{x) 



+ 



b{x) 



sm 



{s,x) 



COS 



ds 



ds + e'R.^t,x) 



u^{x) ■ Vb{x) 



1 — COS 



b{x) 

Finally remembering that due to Lemma ^ 

. f <f>{s,x) \ 
sm ds 

V e y 

with the usual uniform bounds on 7^^, yields Proposition |2 



(s,x) 



ds + en%t,x)^ 









+ 











< en%t,x), 



The asymptotic behaviour of p is obtained in a similar way using the fact that p is propor- 
tionnal to the Jacobian J{t,x) = \ det DX{t, x)\. 

Proposition 3. Consider a function b satisfying assumptions {H0)(H1). Let po be a nonneg- 
ative function in W'-'^'°^{K'^), and uq be a vector-field in iy'''°°(R2) (s > I). For all T < T* 
and £<£t as in TiieoremEl denote by {p, u) the solution of in L°°([0, T],W'-^''^(R'^)) 
and L°°([0,r],VF^'°°(R2)) respectively Then 

p{t, X{t, x)) — pq{x) ^1 + tuQ ■ V log b{x) sin((/)e(t, x)) — Iuq ■ V log b{x) cos((/>e(t, x))^ 
converges strongly to in L°°{[0, T] x R^), where the phase (p^ is defined as previously by i)4.I|) . 
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Proof of Propositions! As long as the solution of (|1.2|) is regular, the equation governing p 
can be rewritten 

— (logp) = V • u, 

where ^ denotes as usual the derivative along the trajectories associated with the flow. Of 
course, the Liouville theorem implies that the equation on the Jacobian of the flow states 

> = ^-"- 

Then, for all e < et, ah t € [0, T] and all x G R^, 

p{t,X{t,x)) = po{x)J(t,x), 

since Jo(t, x) = det(/d) = 1. 



From Lemma m we then deduce that 
(4.2) 

p{t, X{t, x)) — pq{x) ( 1 + tuo ■ V log b{x) sin 



(f){t,x) 



tu^- Vlogbix) 



< Cte. 



Plugging the approximation of the phase obtained previously 

(l){t,x) _ 



(4.3) 



+ e7^^(^,a;) 



back into formula 1)4. 2|) leads then to the expected asymptotics. 

4.2. Inversion of the characteristics. In this section we shall prove Theorems El and ID 
Prom now on T* is the time given by Theorem 121 and we will call T any time smaller than T* 
(in the following we will also suppose e < as given in Theorem IJJ. 

Let X~^{t,x) be the point at time of the trajectory reaching x at time t. By Proposition |2l 
we have 

u{t,x) = uo(X~^(t,x))cos ((j)e(t,X~^{t,x)))-UQ{X''^{t,x))sin (^^(t, X"^(t, x))) +^7^^(^, x) 



with the usual uniform bounds on 7^^. That remainder function Vf{t,x) is liable to change 
from line to line in this paragraph. 

By Proposition n there is a constant Ct (depending on T, uq and 6), such that 
(4.4) Vx€R^ VtG[0,r], \X-'^{t,x)-x\<CTe, 

so we can write rather 

u{t, x) = uq{x) cos {4)£{t, X~'^{t, x))^ — Uq{x) sin (^4)e{t, X^^{t, x))^ + eTZ^{t, x). 

By definition of in 1)4. we have, using again (|4.4() . 

4)e{t,X-^{t,x)) = + ^(X-i(t,x) - x) • V6(x) - tu^{x) ■ Vlog6(x) + e7^^(^,x) 
hence defining 

e,{t, x) =^ x) - x) • V6(x) - tu^{x) ■ V log 6(x) 
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we have 

(4.5) Mt,X'\t,x)) = ^^+6,{t,x)+e'JZ%t,x). 

Now we shall try to make 9^ more precise. According to Lemma Inland the approximation for 
the phase derived in the previous paragraph, we have 

b{x) V 

~^~b(~y ~ [Mt,X~'^{t,x)) +en^{t,x)j^ - etv{x) +e'^n%t,x) 
where again we have used ()4.4() . So we obtain, using the fact that v ■ V6 = 0, 
0e(t,x) = -tno(x) • Vlog6(x)sin (4>e{t,X~^{t,x)) + eTZ^(t,x) 



+tuQ{x) ■ V log 6(2;) cos [J)eit,X''^{t,x)) +e7^^(^,x)j +671" {t,x) 

___ ~ b(x)t 
which by ()4.5p yields directly the result 1)1.4^ . defining 6e{t,x) = 9e{t.,x) H 

Theorem 121 is proved. 

The proof of Theorem 0] is now immediate: we use the formula obtained in Proposition |21 and 
replace pe{t,X{t,x)) by pe{t,x) using the above formulation oi X^^{t,x). The result follows. 
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